Unconventional quasiparticle lifetime in graphite 
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The influence of electron-electron scattering on quasiparticle lifetimes in graphite is calculated. 
In the limit when the Fermi surface is reduced to isolated points in the Brillouin Zone, the suppres- 
sion of screening leads to non Fermi liquid behavior. The inverse lifetime increases linearly with 
energy, in agreement with recent experiments. Similar features should also be present in narrow gap 
semiconductors, and in carbon nanotubes. 

75.10.Jm, 75.10.Lp, 75.30.Ds. 

Recent experiments [jlj show that the inverse lifetime of quasiparticles in graphite increases linearly with energy 
over a broad energy range ~ 0.3 — 4eV. The authors of ref. H argue that this effect may be due to the unusual 
dispersion relation of plasmons in a 3D stack of conducting planes, where the electrons are confined to move in two 
dimensions In the present work, we analyze an alternative explanation of the energy dependence of the lifetime, 
although also associated to the electron-electron interaction. 

Graphite is a semimetal. Band structure || calculations show that electrons are confined to narrow holes near 
the edges of the hexagonal Brillouin zone. The main features around the Fermi level arc well described by a simple 
Hiickel theory Q , which includes only the out of plane tt orbitals located at each carbon atom. Intraplane hopping is 
much larger than the interplane hybridization. If we neglect the coupling between planes, the Fermi surface is reduced 
to two isolated points at the corners of the 2D Brillouin zone. The band dispersion is linear around these points, 
ek = hvp\k\. In units of the hybridization between neighboring tt orbitals, t, we find that Tivf = 3ta/2 where a is 
the distance between carbon atoms. The density of states is zero at the Fermi level, within this approximation. The 
long wavelength properties can be approximated by means of an effective 2D Dirac equation Q , instead of the more 
familiar effective mass approximation, used when the bands are parabolic at low energies. 

Because of the vanishing density of states at the Fermi energy, the Coulomb interaction is not screened in a 
conventional fashion. Combining the band structure described above with the bare Coulomb interaction we obtain 
an effective hamiltonian, for a single plane of graphite: 

n = nv F j M(f)(ia xdx + io ydy )m + ^ J J # ^M*)*ftMV (1) 

where a x and a y are Pauli matrices, and ^ is a two component spinor associated to the two bands which meet at 
the Fermi points. In a honeycomb carbon lattice, t w 2.4 eV. High energy screening processes are included in the 
effective dielectric constant, eo- We expect that ([[J) describes well the physics of a single layer of graphite for energies 
< t. 

The hamiltonian ([!]) defines a renormalizable theory || in the field theoretical sense • Perturbation theory leads 
to logarithmic dependences on the high energy cutoff, which can be incorporated into new, renormalized, parameters 
in a standard way. The existence of logarithmic corrections can be inferred from the fact that the actual coupling 
constant in (Q), e 2 / '{e^fivp) is dimensionless. 

A detailed analysis of the renormalization procedure used to deal with the 2D hamiltonian (Q) is given in ref. [pj. 
That scheme can be generalized to a system of weakly interacting planes. Interplane couplings can be of two types: 

i) Coulomb effects. They increase the number of diagrams that need to be calculated, as an electron in a given 
plane can be scattered by electron-hole pairs in other planes. The low energy properties are, however, not affected. 
The RG scheme proceeds by integrating out slices in energy and momentum space. The momentum transfer in 
interplane processes is bounded by the inverse of the interlayer distance, d -1 . Thus, at sufficiently low scales, the 
electron-electron scattering retains the original dependence on the momentum transfer. The finite corrections induced 
by interplane Coulomb interactions can be included in a manner similar to the scheme used in ref. (pj. The details 
are given below. 

ii) Interplane electron hopping. This interaction is responsible for the 3D features of the bands of graphite. It leads 
to deviations from the linear dispersion relation used to define (Q) at low energies. We expect that these deviations 
will become significant at energies comparable to the interplane hopping, k, 0.27 eV ||. In addition, the density of 
states at the Fermi level becomes finite. Hence, metallic screening takes place at length scales greater than kp^, where 
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kp T = 4ire 2 N(ep). Because of the smallness of N(ep) 0, the associated energy scale, hvpkFT is much smaller than 
the previous one. 

Quasiparticle lifetimes are also influenced by phonons [jlO| , and by low energy, out of plane plasmons [Ell . The 
phonon bandwidth in graphite is ~ 0.20 eV, and the out of plane plasmons have energies ~ 0.05 eV. Thus, for 
quasiparticle energies > 0.2 eV, these processes should give a constant contribution, independent of the quasiparticle 
energy. 

The preceding discussion implies that a straightforward generalization of the hamiltonian ([!]) , including the existence 
of a stack of planes coupled solely by the Coulomb interaction describes well the electronic properties of graphite, in 
the energy range 0.2 eV < e < 2 eV. 

A remarkable feature of the perturbation analysis of hamiltonian ([!]) is that logarithmic divergences appear in the 
corrections to one particle properties, like the self energy, but electron- hole propagators are finite ||. This reflects 
the fact that the divergences are due to the singularity of the interaction, and not to density of states effects. The 
intra BZ edge electron-hole propagator at low energies and momenta is ||: 

f 

Xo(u,q) = (2) 

where xo is purely real for Vp\q\ > to, and purely imaginary otherwise. Thus, electron-hole pairs can only be excited 
if Vp\q\ < lo. This region is shown in fig. (|l|). 

The screened Coulomb potential, including interplane scattering, can be written as 

_ 27re 2 sinh(|g|d) 
*\ w > 1) — i 2 K 6 ) 



cosh(\q)d) + 22jg S mh(\q\d) X o(u>, q) 



e o|g| 

where d is the distance between planes. It is interesting to note that Re v scr has no poles at low energies, unlike for 
the case of a stack of layers with quadratic dispersion, where a plasmon band, lo p oc |g|, was found 0]. We adscribe this 
difference to the fact that, in the present case, electron-hole pairs can only exist for lo > uf|<z|, while in a conventional 
electron gas it is the opposite. 

Using (^|), the inverse of the quasiparticle lifetime can be written as: 

2 f 1 + cos (^fc_jp) / -»\ 

Im£(u;, k) = — ^ Id k Im v scr \ u - e k ,,k - k'J (4) 

where 4>t_c, is the angle between vectors k and k', and we are summing over the two spins. 

The expression (Q) can be interpreted as the probability for a quasiparticle with frequency lo and momentum k to 
decay into a real quasiparticle of energy and momentum k' . Kinematical constraints in the phase space of final 
states imply that ImS(w, k) ^ only if u> < Vf|A;|. This restriction seems incompatible with the phase space available 
for the creation of electron-hole pairs, shown in fig. ([!]), suggesting that there are no channels for quasiparticle decay 
in the model described by (|l|). 

We must, however, consider with care the limit lim w ^ £fc+ Q+ hnE(u>, k). The simplest situation, which can be 
analyzed analytically, is the lowest order diagram shown in figure (Q). For this case, ImSfw, k) drops discontinuously 
to zero at lo — vp\k\. The magnitude of the step is: 

1 / e 2 \ 2 

lim_ t Im£(w,fc) = —(—£— J hv F \k\ (5) 
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The inverse lifetime, defined in this way, increases linearly with the energy of the quasiparticles. The real part of 
the self energy shows a logarithmic dependence on the high energy cutoff needed to define the model (EJ), leading to 
non Fermi liquid behavior 0. The existence of a finite lifetime, despite the kinematical constraints described earlier, 
can be traced back to the divergence of the density of electron- hole pairs in the forward direction, u> — vp\k\, which 
compensates exactly the reduction in the number of states in which the quasiparticle can decay as u> — > t k ■ Setting 
Tivp — 3/2ta, where t — 2.6 eV, a — 1.4 Aand eo — 2.4 [Oj, we find that the constant of proportionality between the 
inverse lifetime and the quasiparticle energy is 0.049 in eV -1 fs _1 . This value compares well to the experimental one, 
0.029 eV^ 1 fs" 1 |Q. Note that the inverse lifetime should be an average of ImE over a finite interval of energies and 
momenta (see below), and expression (^|) is an upper bound to such an average. 
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The behavior of ImE(w, k) as u> — > e&, including the RPA and interplane interactions, can be calculated numerically, 
and it is shown in figure (|j), using expression (||), with d = 3.35 Aand the parameters given above. 

The kinematical constraints discussed earlier arise from the requirement of momentum conservation. In the presence 
of disorder, quasiparticles have a finite spread in momenta. Because of the sharp rise of Im£ away from the line u> = e^, 
this spread leads also to a finite quasiparticle lifetime. This is shown in fig. (^), where Im£(e&, \k\ — Ak) is plotted, 
with Ak = 0.002 A^ 1 . (Afc) -1 corresponds, roughly, to the mean distance beween scattering centers, in our case, 
I ~ 500 A. Other inelastic scatte ring channels, such as phonons, will also contribute to give a spread in momentum 
and energy to the quasiparticles p0| . 

The lifetimes shown in fig. (^) are consistent with the experimental observations Q]. The explanation that we 
propose here also implies non Fermi liquid behavior in other properties of graphite, such as the conductivity or the 
susceptibility. Note, however, that we expect our model to break down at low energies, < 0.2 eV. 

It was argued that low energy plasmons can be responsible for the unconventional quasiparticle lifetimes in graphite 
jij. A layered 2D electron gas has plasmons above a certain threshold which depends linearly on momentum, 
^treshoid = Vpi\q\ f§] • These collective excitations give rise to a new decay channel for quasiparticles with veloc- 
ity v qp = hk/m > v p i. In most cases, decay into low energy plasmons cannot take place at the Fermi level, because 
v p i > vp. It is proposed in ref Q| that, in graphite, v p i < vp. In that case, the low energy plasmons merge into 
the electron-hole continuum. The influence of these hybridized excitations on the lifetimes is rather complicated, and 
does not lead to a precise linear dependence on the energy of the quasiparticles. We believe that, in any case, the 
peculiar band structure of graphite implies that screening is dominated by interband transitions, which preclude the 
existence of in plane low energy plasmons, as deduced from eq.(|J). 

The results depicted in fig.([|) have been obtained by combining a many loop propagator, the RPA bubble modified 
by interplane effects, with a zero loop description of the quasiparticles. We neglect the renormalization of the 
quasiparticle pole, calculated in As we scale towards low energies, the quasiparticle pole loses spectral strength, 
in the manner analyzed in ||. We expect, however, the main result of this paper, shown in fig. (Q), to be weakly 
dependent on this renormalization. Unlike in ID conductors, the hamiltonian (|l|) flows towards a free fixed point [p|, 
which makes plausible the use of unrenormalized quasiparticle propagators in the calculation of the lifetimes. Note, 
however, that the bare coupling constant, e 2 /e^fivp, is of order unity, although most likely reduced by density of 
states factors, ~ (27r) -2 . 

The main physical basis for the unconventional dependence of the quasiparticle lifetimes in energy, and the break- 
down of Fermi liquid theory, lies in the absence of metallic screening. Thus, we expect that other semimetals may 
exhibit similar behavior. In particular, it is well known that the band structure of zero gap semiconductors can be 
approximated by the 3D Dirac equation |l3f| . In these materials, we expect that a description similar to the one used 
here to be valid down to the lowest energies. The graphite structure is the basis of the carbon nanotubes which 
should also exhibit the non Fermi liquid behavior analyzed here. 
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FIG. 1. Region in phase space available for electron-hole excitations. The nature of the interband e-h pairs is sketched in 
the diagram. 



FIG. 2. Lowest order contribution to the quasiparticle lifetime. 



FIG. 3. Imaginary part of the selfenergy as function of frequency for various momenta: Full line, k = 0.1 A 1 , dashed line, 
k = 0.2 A -1 , broken line, k = 0.4 A -1 . 



FIG. 4. Inverse quasiparticle lifetime, as defined in the text, as fucntion of the energy of the quasiparticle. 
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